arXiv:hep-ph/9706245vl 4 Jun 1997 


1 


Beauty hadron lifetime ratios and the problem of t(A;,) 
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We describe a theoretical approach to compute inclusive widths of heavy hadrons, based on a systematic 
expansion in the inverse powers of the heavy quark mass. The method reproduces quite well the experimental ratios 
of B mesons lifetimes. As for T{Ab)/T{Bd), we present a QCD sum rule calculation of the 0(rn^^) contribution 
to T(Ab); we conclude that such correction is not responsible of the small experimental value of r^Ab)/T{Bd)- 


1. INCLUSIVE DECAYS OE HEAVY 
HADRONS 

A simple theoretical way to consider inclusive 
decays of hadrons containing one heavy quark is 
the Spectator Model, which states that only the 
heavy quark participates in the transition, while 
the light degrees of freedom are unaffected by the 
process. As a result, the model predicts that all 
the hadrons with the same heavy quark should 
have the same lifetime. 

On the other hand, the experimental ratios of 
beauty hadron lifetimes are : 

= 1.06 ± 0.04 = 0.99 ± 0.05 (1) 

T\Bd) T\Bd) 

444 = 0.79 ± 0.06 . (2) 

Eq.(||) shows that the ratio T(A;,)/r(Rd) signif¬ 
icantly differs from the Spectator Model predic¬ 
tion. 

One can employ a more refined theoretical ap¬ 
proach, based on QCD, to compute inclusive de¬ 
cay widths of heavy hadrons; the method consists 
in an expansion in powers of ttiq^ |Q. 

By using the optical theorem, the inclusive 
width of a hadron containing the heavy quark Q 
into a final state with assigned quantum numbers 
/ can be written as: 


where T{Q ^ Xf ^ Q) is the transition opera¬ 
tor describing Q with the same momentum in the 
initial and final state: 

f = i J d^x T[CUx)Cl{0)] . (4) 

Cw is the effective weak lagrangian for the decay 
Q ^ Xf. The time ordered product in is com¬ 
puted by an operator product expansion (OPE): 

f = Y, ( 5 ) 

i 

where Oi are local operators ordered by increas¬ 
ing dimension, and Ci are coefficients ordered by 
increasing powers of IjuiQ. As a result, the inclu¬ 
sive width of a heavy hadron Hb (from now on we 
will refer to the beauty sector) is represented by 
a sum of expectation values over Hb of operators 
with increasing dimension: 


T{Hb ^ Xf) = 


GAml 


\KM'(‘ 


+ 


1927r3 
j<Hb\hb\Hb > 

2MHb 

Cg < Hb\bigs(7 ■ Gb\Hb > 
ml 2 Mh, 

< Hb\Of\Hb > , 


E 


2Mm 


o(^ 


\mz 


( 6 ) 


T{Hq ^ Xf) 


2/m < HqITIHq > 
2Mhq 


( 3 ) 


KM is the CKM matrix element relevant in the 
considered decay. 
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The first operator in is bb, which has dimen¬ 
sion D = 3] next, the chromomagnetic operator 
Oc = b^a^j^vG^^b, responsible of the mass split¬ 
ting between hadrons which differ only for the 
heavy quark spin orientation, has D = 5; finally, 
Of have D = 6 and contribute to in 

the expansion. < bb >h^=< Hb\bb\Hb > /^Mh^ 
can be further reduced by using the heavy quark 
equation of motion in the limit rub —*■ oo: 




< Hb\OM > 

2Mh, 


( 10 ) 


fiQ depends on the spin J of Hb'. = —2[J( J-|- 
1) — |]A 2 . It enters with /i^ in the mass formula 
for a heavy hadron in the mb —>■ oo limit: 

+ . ( 11 ) 


< bb >Hh 


< Og >h, 

2mj 


< Otv >Ht 

2ml 



(7) 


Ott = b{iD)'^b is the heavy quark kinetic energy 
operator due to its residual motion inside the 
hadron. 

Some important features should be noticed in 
d)- The first term of the expansion reproduces 
the Spectator Model result; 0(m^^) terms are 
absent [HJJ] since D = 4 operators cannot con¬ 
tribute, being reducible, by equation of motion, 
to the D = 3 operator bb; finally, Oq and Ot^ are 
spectator blind, i.e. they are not sensitive to the 
light flavour. 

terms come from four quark operators, 
which can be classified as follows |^: 


The parameters A, /i^, A 2 are independent of 
the heavy quark mass. For B mesons, can 
be derived experimentally, since it is related to 
the mass splitting: = 3{Mq, — M^)/4 ~ 

0.36 GeV^. On the other hand: ^Q^Ab) = 0, 
since the light degrees of freedom have zero total 
angular momentum relative to the heavy quark. 

Various theoretical determinations exist for 
/i^(i?d). From the mass formula & and using 
experimental data on the mass of the Af, baryon, 
one obtains: 

fj.l{Bd) - A{^b) = 0.002 ± 0.024 GeV^ ( 12 ) 

and hence it can be assumed: ^f^{Bd) — /r^(Ab)[]. 

In order to evaluate the matrix elements of G = 
6 operators, different approaches must be used for 
mesons and baryons. For B mesons, factorization 
approximation gives: 


ces. Their contribution accounts for the presence 
of the spectator quark, which participates in the 
decay through the mechanisms known as weak 
annihilation and Pauli interference (the latter de¬ 
pends on the presence of two identical quarks in 
the final state). 

In order to use the result (^ one has to know 
the matrix elements of all the previously men¬ 
tioned operators. Let us define: 


Mg (77b) — 


< Hb\OG\Hb > 


2M 


Hb 


(9) 


Oy-A — 

bLl^iqLqLltibL 


< Bg\Oy_^\Bg > 

Ol-P = 

bnqLqLbn 


2Mb, 

II 

1 

bLiiJ.t°-qLqLitJ.t‘"bL 


< Bg\Og_p\Bg > 

rpQ _ 

-^S-P ~ 

bRf^qLqLf^bn 

(8) 

2Mb, 

where = 

A“/2, A“ being the Gell-Mann matri- 

< Bg\T^_^\Bg >-- 


f mb + mg \ _ 2 

V Mb, j ~ 


Mf 


fsg being the Bg decay constant. The same tech¬ 
nique cannot be applied to baryons, so that one 
must use other approaches. In the case of Af,, 
heavy quark symmetries reduce the relevant op¬ 
erators to two 01: 


Oy_a = bL-f^bLqLltiQL 
Oy_A = bLl^qLqLl^ibL ■ 


(14) 


^This is confirmed by QCD sum rule estimates 
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One can parameterize the matrix elements as fol¬ 
lows: 

< Oy-A 

< Oy-A 


The parameter b is fixed to 6 = 1 from the 
QCD sum rules analysis of /aj 0, defined by: 
< 0|Jc|Ab(u) >= fAt{'tpv)c ii’v is the spinor for 
a At of four-velocity v). 

The hadronic representation of the correlator 
dil) can be obtained saturating it by baryonic 
states, and considering the double pole contribu¬ 
tion in the variables to and u': 


Jb^b 


48 


= -B 


< Oy j^ >Ai 


2Ma, 


(15) 


in the valence quark approximation: B = 1. 

It should be stressed that only large values of 
the parameter r in (|^) (r = 3 — 4) could explain 
the observed difference between r(Ab) and T{Bd). 
Calculations based on quark model 0 or emloit- 
ing the experimental splitting: 0 give 

values of r smaller than required. In order to have 
a different theoretical estimate, we computed by 
QCD sum rules the matrix elements of the rele¬ 
vant D — 6 operators over At- This calculation 
will be described in the next section. 

2. QCD SUM RULE CALCULATION OF 

< Oy_j^ >At 

The starting point to evaluate < Oy_^ >At by 
QCD sum rules within the Heavy Quark Effective 
Theory (HQET) is the correlation function: 


UcoioJjOj') = (1 + 

<0|T[Jc(x)0^_^(0)Jr5(2/)]|0> (16) 

where C, D are Dirac indices; the variable w {lo') 
is related to the residual momentum of the incom¬ 
ing (outgoing) baryonic current = mtv^ + 
with = ujv^. A suitable interpolating field for 
At, in the limit mt —> oo, is given by 0: 

Jc{x) = e"^'‘{q^\x)TTq^{x)){h^)cix) (17) 

where T means the transpose, i, j, k are colour 
indices, hy is the effective heavy quark field, 
r = ( 775(1 + b'^) and C is the charge conjuga¬ 
tion operator, r is the At light flavour matrix: 




(Aa6 - w)(Aa6 - w') 


-f ... 


(19) 


at the value uj = uj' = Aa,,. 


Aa^ represents the 


binding energy of the light degrees of freedom in 
At'. Ma^ = rrit + Aa^ and must be derived within 
the same QCD sum rule framework. On the other 
hand, in the Euclidean region, for negative values 
of UJ, uj', the correlator (^) can be computed in 
QCD, in terms of a perturbative contribution and 
of vacuum condensates. The result can be written 
in a dispersive form: 


n 


OPE 


{uj,uj') = J da da' 


Pn(o-, cr') 


{a — uj){a' — uj') 


( 20 ) 


where possible subtraction terms have been omit¬ 
ted. The explicit expression of the computed 
spectral function pn can be found in |0|. 

A sum rule for < Oy_j^ > can be derived by 
equating and H'^'^^. Moreover, invoking 

a global duality ansatz, the contribution of the 
higher resonances and of the continuum in 
in ( p^ can be modeled as the QCD term outside 
the region 0 < a; < ujc, 0 < w' < uJc, with uJc an 
effective threshold. Einally, the application of a 
double Borel transform to both H*^^^ and 
in the momenta uj,uj': 


B{E,) 


1 


a — UJ 


B{E^) 



E, 

1 


( 21 ) 


r = 


1 

71 



(18) 


(and similar for uj' with the Borel parameter E 2 ) 
allows us to remove the subtraction terms in (^^ , 
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Figure 1. Sum rule for < Oy_j^ >A|, as a function 
of the Borel variable E. 


that are polynomials in oj or uj' (the Borel trans¬ 
form of a polynomial vanishes). Moreover, the 
convergence of the OPE is factorially improved 
by the transform, and the contribution of the 
low-lying resonances in is enhanced for low 
values of the Borel variables. The symmetry of 
the spectral function in ct, a' suggests the choice 
El = E 2 = 2E. The final sum rule reads: 


/■““ /■““ cr -I- cr' 

^ Jo Jo - 

The result of the numerical analysis of (|2|) is 
depicted in Fig. |l|. A stability window is ob¬ 
served, starting at E ~ 0.2 GeV and continuing 
towards large values of E; in this range the re¬ 
sult for {Oy_^)A^ is independent of the external 
parameter E. The variation of the result with E 
and with ojc provides us with an estimate of the 
accuracy of the numerical outcome. We find [Tl| : 

{01_a)a, ^ (0.4 - 1.20) X 10-3 GeV^ , (23) 

a result corresponding to r ~ 0.1 — 0.3. 

As for B in (|^, since our computational 
scheme considers only valence quark processes, 
a sum rule for the matrix element in ( p^ would 
give B = 1. The explicit calculation confirms this 
conclusion. Using the computed r and B in the 
formulae in ref. we get: T{Ab)/T{Bd) > 0.94. 


3. CONCLUSIONS 

We presented a method based on the 1 /toq 
expansion to compute inclusive widths of heavy 
hadrons. We also described the QCD sum rule 
calculation of the matrix element < Oy_j^ >Ab 
contributing to 0{m^^) to the A;, lifetime. The 
result: T{Ab)/T{Bd) > 0.94 implies that such cor¬ 
rection cannot explain the observed difference be¬ 
tween and T{Bd). If the experimental data 

will be confirmed, a theoretical reanalysis of the 
problem would be required, as suggested for ex¬ 
ample in ||l^ , |l^ . 
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